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The 7.6 eV electromagnetic transition between the nearly degenerate ground state and first excited 
state in the 229 Th nucleus may be very sensitive to potential changes in the fine-structure constant, 
a — e 2 /he. However, the sensitivity is not known, and nuclear calculations are currently unable to 
determine it. We propose measurements of the differences of atomic transition frequencies between 
thorium atoms (or ions) with the nucleus in the ground state and in the first excited (isomeric) 
state. This will enable extraction of the change in nuclear charge radius and electric quadrupole 
moment between the isomers, and hence the a-dependence of the isomeric transition frequency with 
reasonable accuracy. 
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I. INTRODUCTION 

The isotope 229 Th has the lowest known excited state 
in nuclei; recent measurements show that the 3/2+ state 
lies just 7.6 eV above the 5/2+ ground state The 
width of this level is estimated to be about 1CP 4 Hz [2] 
which may explain why it is so hard to find the direct 
radiation in this very weak transition. Nevertheless, the 
frequency is within the range of lasers, and it has been 
proposed to use this narrow nuclear transition as possible 
reference for an optical clock of very high accuracy 
Additionally, this transition could be a sensitive probe of 
possible variation of fundamental constants [3| because 
the near degeneracy of these isomers is a result of can- 
celation between very large energy contributions (order 
of MeV). Since these contributions would have differ- 
ent dependences on fundamental constants, any varia- 
tion would be enhanced in the transition frequency. In 
Ref. [i]], the relative effects of variation of a and the 
dimensionless strong interaction parameter, m q / Aqcd 
were estimated to be enhanced by 5 orders of magnitude. 

An enhancement to a-sensitivity of this magnitude 
would have very important consequences for laboratory 
searches of a-variation. Because the isomeric 229 Th reso- 
nance has a narrow linewidth and an extraordinary insen- 
sitivity to external perturbations, an optical clock utilis- 
ing this reference may have very high accuracy and high 
immunity from systematic frequency shifts \3Q. By com- 
paring this "nuclear clock" frequency with that of any 
other narrow optical or microwave transition (e.g. the 
Cs or Hg + frequency standards) one can test variation of 
fundamental constants. Coupled with the enhancement 
in sensitivity, such a set up would be the most sensitive 
laboratory probe of a- variation to date, possibly gaining 
several orders-of-magnitude improvement over the cur- 
rent limits of a/a = (-1.6 ± 2.3) x 1CT 17 year" 1 |[. 

The sensitivity of the transition frequency to variation 



of a can be expressed as 
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where AVc is the difference in Coulomb energies be- 
tween the two isomers, and K is the enhancement factor: 
K = AVc/u. Since the Coulomb energy of this nucleus 
is of order 10 9 eV, even a relatively small variation in Vc 
could produce a large enhancement. For AVc — 100 keV 
and Sa/a = HT 16 , Sui = HT 11 eV = 2.4 x 10 3 Hz which 
is 4 orders of magnitude larger than the limits placed on 
shifts in atomic transitions in Ref. [5[. 

However, different nuclear calculations give wildly dif- 
ferent values for AVc- Refs. [1,0] claim that both iso- 
mers have identical deformations and therefore the same 
Coulomb energies to within roughly 30 keV (correspond- 
ing to K < 4000). Ref. jl gives a value of 30 keV, 
while the calculations of [9( give values in the range 
-300 keV < AVc < 450 keV, depending on particulars 
of the model used. Lastly, Ref. [Ic| uses Nilsson wave 
functions to show that the value of AVc a s a function of 
deformation changes from 1.5 MeV at zero deformation 
down to -0.5 MeV at 8 = 0.3. They conclude that a very 
small value of the Coulomb energy shift seems improba- 
ble. 

In this paper we propose a different method for extract- 
ing sensitivity to a- variation using direct laboratory mea- 
surements of the change in nuclear mean-square charge 
radius, A(r 2 }, and electric quadrupole moment, AQo, 
between the isomer and the ground state nucleus. In 
Section|H]we present a simple geometric model of the nu- 
cleus to relate the observable nuclear parameters to AVc 
and hence K . We show that this model is self-consistent 
by comparing to the nuclear calculations of Ref. 0] . 

Once the change in Coulomb energy has been mea- 
sured, the change in nuclear energy AEn will be known 
also, since they almost cancel for this transition: AEn + 
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AVc = 7.6 eV. This change in nuclear energy can be in- 
terpreted in terms of variation of the dimensionless ratio 
m q /h.QCD where m q is the light quark mass and K-qcd is 
the pole in the running strong coupling constant. Vari- 
ation of this ratio will also be enhanced in the 229 Th 
transition Q ; calculations may be found in [1, [ll| . 

The change in mean-square nuclear radius can be ex- 
tracted using the isomeric field shift for an atomic tran- 
sition. In principle, any transition in any 229 Th ion or 
the neutral atom can be used. There are two approaches. 
The first is entirely empirical: by combining the measure- 
ments of isomeric shifts and isotopic shifts for the same 
transition, one can extract the ratio of A(r 2 } for the iso- 
mer to the isotopic change in mean-square radius. The 
second approach does not require the additional measure- 
ment of isotope shift, but it does require high-precision 
atomic calculations. We provide more details and neces- 
sary calculations in Section InTl 

To extract the change in nuclear quadrupole moment, 
the hyperfine structure may be used. The hypcrfinc 
structure constant B, which can be determined exper- 
imentally, is proportional to the quadrupole moment Q. 
Therefore one must measure B for both the ground state 
and isomeric 229 Th. The value of Q for the ground state 
is known to within 20%; better accuracy can be obtained 
using the calculations presented in Section llVl 

The radiative lifetime of the metastable 229m Th nu- 
cleus is estimated to be a few hours Q, however this 
may be reduced if the energy of the excited state exceeds 
the ionization potential since an electron autoionization 
channel may open. The successive ionization energies of 
thorium ions are 0: 6.3 eV (Th I), 11.9 eV (Th II), 
20.0 eV (Th III), 28.8 eV (Th IV). Therefore, the atomic 
experiments are likely to be easier for ionized thorium 
since the ionization energies exceed the excitation energy. 
In fact, Th IV may be the best choice since it is alkali- 
like and is amenable to laser cooling and trapping 
This ion has the additional advantage that calculations 
are likely to be more accurate, although we stress that 
we can obtain reasonable accuracy with any ion that ex- 
perimentalists may find convenient. 



II. 



GEOMETRICAL NUCLEAR MODEL 



In this section we use a simple geometric model to re- 
late the Coulomb energy of a nucleus to the experimen- 
tally observable mean-square charge radius and electric 
quadrupole moment. We assume that both the ground- 
state nucleus and the lowest-energy isomer are uniform, 
hard-edged, prolate ellipsoids. Let a and c be the semi- 
minor and semi-major axes, respectively. We define Rq 
as the equal-volume spherical radius and £ as the eccen- 
tricity, so that 



Rf. = a 2 c and ( 2 
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We can extract from experiment the mean-square ra- 



dius and intrinsic electric quadrupole moment (see Sec- 
tions [ni] and QV} defined as 

(r 2 ) = / r 2 p(r)d 3 r (3) 

r 2 (3 cos 2 (6>) - l) p(r)d 3 r (4) 



where p(r) is the electric charge density normalised to 
unity. The intrinsic quadrupole moment is related to the 
laboratory quadrupole moment of the ground rotational 
mode by (see, e.g. fl3j ) 
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'(7+l)(27 + 3) ' 
For our hard-shell prolate nucleus, one finds 

(r 2 
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- (2a 2 + c 2 ) and Q Q = - (c 2 - a 2 ) 
5 5 



We wish to express the Coulomb energy in terms of 
these measurable quantities. Using formulas presented 



14| we find 
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and in terms of (r 2 )and Qo we finally obtain 
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With this equation we can extract AVc if we know A(r 2 ) 
and AQo between the 229 Th isomers. Note that Vc and 
AVc a re vastly more sensitive to changes in (r 2 ) than 

Qo- 

To estimate the effect of skin thickness, we use a spher- 
ical Fermi distribution model: 



p( r ) = Po 1 + exp 



r-C 



The Coulomb energy is 



V c = 
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and one sees that Vc is not sensitive to the skin thickness 
parameter z. 

With the current data for 229 Th, r rms = 5.6807 ± 
0.0509 fm [H and Qi ab = 4.3 ±0.9 eb 0, we obtain 
Qo = 13.4, an eccentricity £ 2 = 0.440, and Coulomb 
energy Vc = 967 MeV. We estimate z = 0.5 fm. In this 
case the change in Coulomb energy can be expressed 

AV C A(r 2 ) AQ Az , s 

— A- = -506 + 23 + 17 — 10 

MeV) (r 2 Qo z 
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TABLE I: The values of AQo, and V c 

are reproduced from Ref. [y] and used to calculate the value 
of AVa shown in the last line using our simple geometrical 
model. SkM* and SIII refer to two different energy function- 
al, while HF and HFB refer to Hartree-Fock and Hartree- 
Fock-Bogoliubov, the latter includes pairing correlations; for 
details see Ref. [Si. 



SkM* 
HF HFB 



SHI 
HF HFB 



r rms (fin)" 
Qo (fm 2 ) a 
Ar rms (fm) a 
AQo (fm 2 )° 
V c (MeV) 6 

AV C (MeV) 6 
AVc (MeV) 



5.7180 5.7078 
9.5461 9.3717 
-0.0038 0.0039 
-0.1824 0.2756 
924 925 



5.7817 
9.3542 
0.0000 
-0.0339 
912 



5.7769 
9.1643 
-0.0005 
-0.0495 
912 



0.451 
0.419 



-0.307 
-0.327 



-0.098 0.001 
-0.036 0.029 



"From Ref. ||, Table II. 
6 From Ref. Q, Table I. 



from which the sensitivity of the transition to a- variation 
is easily deduced. Note that the contribution of Az is 
small. 

As a consistency check, we have recalculated AVc us- 
ing the values of r rms , Ar rms , Qo, and AQo calculated 
in Ref. [9(; this is shown in Tabled] That we are able to 
reproduce their results shows the validity of geometrical 
model. The differences seen in the SIII entries of Table U 
(last two columns) are probably due to lack of numer- 
ical precision. If measurements of Ar rms and AQo are 
made accurately, the model should suffice even when the 
measurable nuclear parameters are small. 



III. MEAN-SQUARE RADIUS 

In the previous section we showed that A(r 2 ) is the 
most important quantity for determining AVc and hence 
sensitivity to a-variation. In this section we show how 
A(r 2 )can be extracted from the isomeric shift of any 
atomic transition, obtained by comparing 229 Th and 
229m Th. This is similar to the usual isotope shift. 

The shift in energy of any transition in an isotope with 
mass number A' with respect to an isotope with mass 
number A can be expressed as 



Av A '' A = (fc N MS 



fcsMs) FA ( r ) 



„2\A',A 



(11) 

Here the first term is the "mass shift" due to the finite 
mass of the nucleus and the second term is the "volume" 
or "field" shift due to the finite size of the nuclear charge 
distribution (see, e.g. In the case of the isomeric 

shift that we are interested in, the mass shift vanishes 
since isomers have equal mass. Thus in order to extract 
A(r 2 ) from a measurement of isomeric shift Av m for an 
atomic transition we need simply divide by the field-shift 



TABLE II: Calculated energies and field shift constants of 
transitions in Th IV. The last column shows expected "order 
of magnitude" isomeric shifts in 229 Th, assuming |Ar rms | = 
0.004 fm, however the actual shift could differ by an order of 
magnitude. All transitions are to the 5/5/2 ground state. 



uj (cm" 1 ) F \S^\ 

Level Expt. Calc (GHz/fm 2 ) (GHz) 

0.09 
1.4 
1.5 
6.3 
2.5 
2.1 



2 (2) 
33 (8) 
35 (8) 
146 (4) 

7pi /2 60239 63051 57 (3) 
7p 3/2 73056 76319 49 (2) 



5/7/2 4325 4899 

6d 3/2 9193 11721 

6d 5 /2 14486 17534 

7s 1/2 23131 24740 



constant F: 



Av m =FA(r 2 ) 



(12) 



These may be calculated or extracted from known isotope 
shifts. 

In Tables lHllIIH and llVl we present calculated field shift 
constants for transitions in several ions of Th. In Table HTl 
we have included an estimated size of the isomeric shift, 
Sv m , assuming that Ar rms = 0.004 fm, which is the mag- 
nitude of the largest shifts in (from the SkM* nuclear 
energy functionals). 

We calculate the field shift constants F using methods 
developed in previous works [13] • Briefly, we perform an 
energy calculation several times, modifying the nuclear 
radius in our codes. F is extracted from the gradient: 



F = dE/d{r 2 ) at r r 



5.6807 fm. 



Calculations of the energies are slightly different for 
a single-valence-electron ion (Th IV) and for two- and 
three- valence-electron ions (Th III and Th II). In the for- 
mer case we use the correlation-potential method devel- 
oped in Ref. [l8| . The second-order correlation correction 
potential E^ 2 ) is used to calculate Brueckner orbitals for 
the states of the valence electron. This techniques takes 
into account dominating relativistic and correlation ef- 
fects and leads to good agreement between theoretical 
and experimental energies as illustrated in Table HT1 

For ions with two and three valence electrons we use 
the combination of the many-body perturbation theory 
and the configuration interaction technique (CI+MBPT, 
Refs. [l9|, [20] ) - The same single-electron correlation 
correction operator Si is used for all three ions, in- 
cluding the single-electron ion Th IV. However, an ex- 
tra two-electron correlation correction operator S2 is 
needed for ions with more than one valence electron (see 
Refs. [Hill for details). The accuracy of these calcula- 
tions is also high, as is illustrated in Tables HTT1 and HVl 

For Th II there are experimental isotope shifts avail- 
able [2l| and we compare them with our calculations in 
Table IIVI Note that the mass shift has been ignored 
here: while fcsMS is difficult to evaluate accurately, /cnms 
is easily extracted from the transition frequency and is 
proportional to it. If we assume that &nms and fcsMS 
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TABLE III: Calculated energies, w (cm^ 1 ), field 
shift constants, F (GHz/fm 2 ), and isotope shifts, 
^232,230 (lo-Scm" 1 ), of some transitions in Th III. All 
transitions are to the 5/6d 3 H4 ground state. Note that, 
while we believe the 6d 2 3 F 3 , 6d 2 3 F 4 , and 6d7s 3 D 3 
transitions are accurate, the others are estimates only. 



Level 




to (cm i ) 


F 


^232, 23U 


Term 


J 


Exp. 


Calc 




Calc. a 


6d 2 3 F 


3 


4056 


4023 


24 


165 


6d 2 3 F 


4 


6538 


6795 


22 


147 


&d7s 3 D 


3 


9954 


9204 


118 


804 


6d 2 X G 


4 


10543 


11051 


8 


56 


5/ 2 3 H 


4 


15149 


13358 


-11 


-77 


5/ 2 3 H 


5 


17887 


16068 


-20 


-136 


5/ 2 3 F 


3 


20840 


19080 


-18 


-122 


5/ 2 3 F 


4 


21784 


20366 


-15 


-101 


5/ 2 X G 


4 


25972 


25269 


10 


-66 


5/7P (f , |) 
5/7P (|,|) 


3 


33562 


33402 


13 


92 


3 


38432 


38617 


15 


101 



A(r 2 ) = 0.205 fm 2 , from Ref. [TJ] 



TABLE IV: Calculated energies, to (cm^ 1 ), field 
shift constants, F (GHz/fm 2 ), and isotope shifts, 
^232,230 ( 10 -3 cm -i) ) f som e transitions in Th II. All 
transitions are to the 6d 2 7s J = 3/2 ground state. 

Llvel to (cm' 1 ) F <5!/ 232 ' 23U 

Configuration J Exp. Calc Exp. Calc." Calc. 6 



5/7s 2 


2po 


5/2 4490 4856 


4 


54 


47 


43 


5/6d7s 


4 F° 


3/2 6691 7487 


-53 


-362 


-401 


-362 


5/6d7s 




5/2 7331 8325 


-53 


-365 


-405 


-365 


5/6d7s 


4 G° 


5/2 9585 10045 


-55 


-375 


-406 


-366 


5/6d7s 


4 H o 


5/2 10673 12168 


-53 


-361 


-406 


-367 


5/6d7s 


2 D o 


3/2 11576 13054 


-54 


-367 


-408 


-368 


5/6d7s 


4 D o 


1/2 11725 12897 


-67 


-456 


-460 


-415 


5/6d7s 


2po 


5/2 12472 14564 


-58 


-399 


-463 


-418 


5/6d7s 


4po 


3/2 12902 14233 


-58 


-395 


-444 


-400 


5/6d7s 


4 G° 


1/2 14102 15853 


-79 


-539 


-610 


-550 



a A(r 2 ) = 0.205 fm 2 , from Ref. [TJ 
b A{r 2 ) = 0.185 fm 2 , best fit value. 



are of the same order, then A;nms(1/232 — 1/230) «2x 
10 -8 v is negligible. The second-last column of Table ITVl 
is a calculation with A( r 2)232,230 _ q.205 (30) fm 2 [nf7 
The last column gives values of the isotope shift with 
A(r 2 ) 232 - 230 = 0.185: this is the value that gives the best 
fit of our calculated isotope shifts to the experimental 
data. 

The field shift constant is generally larger for transi- 
tions involving a change in the s-wave configuration, e.g. 
5/5/2 ~ > 7si/2 transition in Th IV and the 5/6d 3 H4 — > 
6d7s 3 D3 transition in Th III. Measurement of the iso- 
meric shift may be easier for these cases. However if 
there are good reasons to use transitions with smaller 
shifts (e.g. the higher-energy transitions in Th III), then 
we recommend the experimentalists contact us for more 
precise values of the constants. Again we stress that 



TABLE V: Calculated electric-quadrupole hyperfino- 
structure constants B for some low energy states of Th IV. 
In the last column, the nuclear electric quadrupole moment 
Q is taken to be 4.3 b. 



Level 


B (MHz) 


5/5/2 


740 Q 


3180 


5/7/2 


860 Q 


3700 


6^3/2 


690 Q 


2970 


6^5/2 


860 Q 


3700 


7P3/2 


1810 Q 


7790 



these constants may be extracted from measured isotope 
shifts with accuracy limited by knowledge of the isotopic 
change in mean-square radius, A(r 2 ) A ' . 



IV. ELECTRIC QUADRUPOLE MOMENT 

Although we have shown in Section [ill that the change 
in Coulomb energy of the 7.6 eV transition in the 229 Th 
nucleus is far more sensitive to (r 2 ) than Qo, AQo could 
still be important if A(r 2 ) is found to be very small. For- 
tunately AQo can be extracted from measurements of the 
hyperfine structure of the isomers by using ([5]) and not- 
ing that Qi a b is proportional to the electric-quadrupole 
hyperfme-structure constants B. 

Since the electric quadrupole moment of the ground 
state 22 Th nucleus is known to about 20% accuracy 
(Qiab = 4.3(9) eb [Hj]), AQo can be extracted by mea- 
suring the electric quadrupole hyperfine structure of both 
isomers. This can be done for any states of any tho- 
rium ion or neutral atom and no atomic calculations are 
needed for the interpretation of the results. 

If better than 20% accuracy is required, the values of 
Qiab can be found by comparision of the calculated and 
measured B. Calculations with this level of accuracy 
for many- valence-electrons are difficult, but can be per- 
formed if required. In this work we present the calcu- 
lations of B for the single-valence-electron ion Th IV. 
The calculations are done with the correlation potential 
method which takes into account dominating correlation 
corrections [18J . The constant B for a particular valence 
state v is found as a matrix element 

B v = A(tf r \\F + SV\\4>?), (13) 

where A is a numerical constant, -0^ r is the Brueckner 
orbital for the valence state v, F is the operator of the 
nuclear electric quadrupole moment and SV is the cor- 
rection to the atomic self-consistent potential due to the 
effect of nuclear quadrupole electric field on atomic elec- 
trons. The same Brueckner orbitals are used as in the 
previous section. The results are presented in Table |V| 
accuracy is expected to be at the level of a few per cent. 
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V. CONCLUSION 

We have presented a simple geometrical model which 
allows one to calculate changes in the Coulomb energy be- 
tween the different isomers given small changes in mean- 
square radius and quadrupole moment; with current data 
the change is given by Eq. (jXOJ) . These parameters can 
be obtained by measurement of the atomic spectra of 
229 Th and its isomer. From the change in Coulomb en- 
ergy, the sensitivity of the isomeric transition frequency 
to a-variation can easily be deduced. 

Two approaches have been proposed for measuring the 
change in mean-square charge radius: in the first the iso- 
tope shift must be measured in conjunction with the iso- 



meric shift. In the second approach measurement of an 
isotope shift is not needed, but atomic calculations are re- 
quired to interpret these measurements. We have shown 
that we can calculate the relevant parameters: namely F 
for extracting A(r 2 } (Eq. [T2]) and B for extracting AQo 
(Eq. [T3|) . We recommend that experimentalists contact 
us for more accurate calculations for the atomic transi- 
tions that they intend to exploit. 
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